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ABSTRACT 

Tidally distorted accretion disks in binary star systems are subject to a local hy- 
drodynamic instability which excites m = 1 internal waves. This instability is three 
dimensional and approximately incompressible. We study the global aspects of this 
local instability using equations derived under the shearing sheet approximation, where 
the effects of the azimuthal variation along distorted orbital trajectories are included 
in source terms which oscillate with local orbital phase. Linear analyses show that the 
excitation of the instability is essentially local, i.e. insensitive to radial boundary con- 
ditions. The region of rapid growth feeds waves into the region of slow or negligible 
growth, allowing the instability to become global. The global growth rate depends the 
maximum local growth rate, the size of the rapid growth region, and the local group 
velocity. We present an empirical expression for the global growth rate. We note that 
the local nature of the instability allows the excitation of waves with m ^ 1 when the 
local growth rate is large. 
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1. Introduction 

The evolution of accretion disks is governed by the transport of angular momentum outward, 
which allows mass to flow inward (see Pringle f98f and Frank, King, & Raine f992 for review and 
further references). In the standard a-disk model (Shakura & Sunyaev 1973; Lynden-Bell & Pringle 
1974), it is assumed that angular momentum transport is due to a local parameterized viscosity of 
the form 



where H and i7 are the disk thickness and the orbital frequency respectively and a is a constant 
of order unity. Since microscopic molecular and radiative viscosities are normally very small, 
turbulence is often invoked as the source of the local viscosity. Consequently, local instabilities 
capable of inducing turbulence play a critical role in the theory of accretion disks. 

Recently, Goodman (1993, henceforth G93) identified a hydrodynamic instability caused by 
tidal distortions of accretion disks. This instability follows from earlier work showing that laminar 
fiows with high symmetry following rectilinear and circular paths may be stable against the growth 
of linear perturbations, but small departure from symmetry may lead to instabilities (Bayly, Orszag, 
& Herbert 1988 and references therein). The original work on this idea was done in the context of 
inviscid incompressible elliptical vortices, so the instability is called the elliptical instability (Gledzer 
et al. 1974; Gledzer et al. 1975; Pierrehumbert 1986; Bayly 1986; Craik & Criminale 1986; for 
further references, see Lubow, Pringle, & Kerswell (1993)). Since the angular momentum of the 
excited waves has the opposite sign of the specific angular momentum of the gas in the disk, and 
of the stars in the binary system, this process leads to the transfer of angular momentum from the 
disk to the binary system. In that role it is important in limiting the outer radius of the accretion 
disk and in providing a sink for the angular momentum transported outward from small radii. 

At smaller radii within the disk the direct excitation of waves by the elliptical instability 
becomes weaker, and eventually completely ineffective. Nevertheless it may still play an important 
role in angular momentum transport in the inner parts of the disk. In unmagnetized disks the inward 
traveling waves excited at large radii can redistribute angular momentum, driving the evolution of 
the disk (Vishniac & Diamond 1989). In magnetized disks the local transport of angular momentum 
is probably mediated by a general magnetohydrodynamic instability first discovered by Velikhov 
(1959) (see also Chandrasekhar 1961) whose importance in the context of accretion disks was 
pointed out by Balbus & Hawley (1991). This process requires only that the disk be magnetized, 
with a magnetic pressure not greater than the ambient gas pressure, and with an angular velocity 
that decreases away from rotation axis. Balbus and Hawley proposed that this process is capable of 
driving a dynamo with a growth rate comparable to i7, a conjecture supported by recent numerical 
calculations (Brandenburg et al. 1995; Hawley, Gammie, & Balbus 1995). One of us (ETV) 
has argued elsewhere (e.g. Vishniac & Zhang (1995)) that this conclusion is inconsistent with 
the evidence for variations of a with temperature and radius found in phenomenological models 
of dwarf novae and X-ray transients. In the absence of a local dynamo, the waves excited by the 
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elliptical instability are also capable of driving a dynamo (Vishniac, Jin, & Diamond 1990, Vishniac 
& Diamond 1992). 

The linear analysis of the elliptical instability in accretion disks was originally performed by 
G93. In binary systems, disks are tidally distorted in a predominantly m = 2 pattern. These tidal 
distortions can destabilize m = 1 inertial waves, whose restoring force is due to the gradient of 
the specific angular momentum in the disk. In real accretion disks the relevant modes are strongly 
affected by the vertical entropy gradient as well and are frequently referred to as internal waves, by 
analogy with underwater ocean waves. G93 found that the local linear growth rate at small radii 
(r < rb) is 

15 Ms 9l 



4 Mp + Ms 1) ' 



where Mp and Mg are the mass of primary and secondary stars, i7b = [G{M^ + M^) I r^^'^ is the 
orbital frequency, and rb is the binary semi-major axis. This formula can underestimate the growth 
rate when r/rb is not small. 

The nonlinear behavior of the instability was studied numerically by Ryu & Goodman (1994, 
henceforth RG) using a quasi-axisymmetric approximation in the meridional {rz) plane. For pertur- 
bation wavelengths small compared with the disk thickness, unstable modes were shown to remain 
nearly incompressible even when nonlinear. The modes saturated chaotically at perturbation ve- 
locity 

^^^(11 ±2)^, (3) 

where 27r/A; is the largest unstable wavelength in the simulations. The two-dimensional power 
spectrum of velocities after saturation was roughly isotropic and extended over a broad range of 
scales in an approximately power-law fashion. 

The significance of this instability hinges on its persistence in the face of realistic damping 
mechanisms. This is dependent both on the local growth rate, and on the sensitivity of the 
growth rate to boundary conditions. Lubow, Pringle, & Kerswell (1993) argued that the elliptical 
instability is likely to prove insignificant on both grounds. They calculated the growth rate for large 
vertical wavelength modes, including the effects of vertical structure and stratification, and derived 
growth rates a factor of ten smaller than the asymptotic formula given in equation (2). They also 
estimated the effects of energy losses at disk boundaries, and concluded that these could result in 
the suppression of the instability. Subsequent work by Vishniac & Zhang (1995) has shown that 
the inclusion of realistic vertical structure does not substantially reduce the growth rates, and that 
the reduction found by Lubow et al. actually followed from their neglect of radial gradients in the 
streamline ellipticity. 

In this paper, we study the global aspects of the local elliptical instability in accretion disks 
using a linear analysis. In §2, we review the linearized equations used and introduce models and 
assumptions. In §3, we describe the results of the linear analysis focusing on global properties of 
the instability and the role of radial boundaries. Summary and conclusions follow in §4. 
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2. Equations and Methods 

We are concerned with small perturbations of a basic state that is stationary in a frame 
corotating with the binary and centered on the primary star. The streamlines of the basic flow are 
closed in this frame, but they are noncircular because of the tidal field exerted by the secondary 
star. Following G93 and RG, we adopt a generalized form of the shearing sheet approximation 
(Goldreich & Lynden-Bell 1965; Julian & Toomre 1966). The fiuid equations are written in a 
pseudo-Cartesian coordinate system xyz whose origin orbits with the basic fiow along a fiducial 
streamline. Spatial gradients of the basic fiow such as shear and vorticity are retained to the lowest 
nontrivial order in the size of the patch divided by the mean radius of the fiducial streamline (tq). 
Strains associated with the ellipticity of the streamlines are included to first order in the tidal field. 

It may be asked how we can adopt a local approximation such as the shearing sheet when 
we are concerned with global questions such as radial boundaries and radial variations in the tidal 
strain. The answer is that we do not want to model any particular disk in realistic detail, but 
rather to explore the effects of generic radial inhomogeneities on the growth rates. This we can 
accomplish by imposing artificial but rapid radial variations in the tidal terms of the shearing sheet 
equations. 

Absent such artificial variations, the fiuid equations are (RG) 

Dp ^ f dvx , dvy ^ dv 



z), ■ " life + 97 + arj = 

1 dp ( 1 dp 1 dp\ 

-fZT + -T, = -2fl {c^^v^ + c^yVy) + b^^-— + b^y-— (5) 

Dt pox \ pox poyj 

DVy l^P , r>/r> , r>r> / , \ , f I. ^ , I. ^ 



JJ^ ' pQy+'^i^ + ^)'^X= -'i^{Cy^V^ + CyyVy)+ [by^-—+byy-—] (Q) 



Dt 



~^ ^ = -^^iCzzVz + bzz TT- (7) 

poz oz \poz oz J 



Dp ^ ( , dvy ^ 



where 



D d d d d 

Dt Ot ox oy Oz . 



Here i7 is the average angular velocity of the basic fiow as it travels along the fiducial streamline, 

and 

is the average shear. 

The coefficients Cij and bi, and the the enthalpy of the basic state, 

Wo{t) = Mt) + f2{t)z\ (11) 
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describe the local effects of the tide on the streamlines. In the absence of the tide, the c's and 
5's would vanish, and wq = (i7^/2)(i7^ — z'^). The time t upon which these coefficients depend 
is really orbital phase along the fiducial streamline. With a predominantly m = 2 tidal pattern, 
the coefficients (cij(t), bi(t), fo(t), and /2(0) oscillate sinusoidally at frequency 2(i7 — i^b), with an 
amplitude that depends on (Mg/Mp) and (ro/rb). 

Thus t plays a double role in our shearing sheet equations. For the basic state it is just orbital 
phase. But for perturbations of that state it is a true time variable, since such perturbations need 
not be periodic. This duality is consistent with our local approximation. 

When linearizing equations (4) to (9), we make a couple of simplifications. 

First, all derivatives with respect to the azimuthal coordinate (y) are dropped, but all three 
components of velocity are included. In the linear regime, as explained in §§1&4, the m = 2 tidal 
pattern excites primarily m = 1 internal waves. The azimuthal wavelength of these waves (27rro) 
is effectively infinite compared to the size of our local coordinate patch. 

Second, the vertical structure of the stationary disk has been neglected by dropping the terms 
involving dwo/dz. Local analysis (G93, Vishniac & Zhang (1995)) indicates that these terms are 
not important in determining the linear growth rate except the radius at which 

(n-n^f = {^ + l)n\ (12) 

where a vertical pulsation mode resonates directly with the tidal force (Lubow 1981). On the other 
hand, this simplification makes the mathematical analysis much easier by permitting a Fourier 
decomposition along the vertical (z) direction. Vertical structure and vertical boundary conditions 
have been considered previously by Lubow, Pringle, & Kerswell (1993) and Vishniac & Zhang 
(1995). Our focus here is on radial boundary conditions and radial variations in the tidal terms. 
We therefore assume perturbations of the form 

Sq(x,z,t) = Sq(x,t)ej''-'. (13) 

Since we want to study radial inhomogeneities in the basic state, we do not assume a Fourier 
decomposition along x. 



Our linearized perturbation equations are 

1 dSp dSvx 



Po dt dx 



ikjvz = 0, (14) 



2^lSvy = -21) {c^^Sv^ + c^ySvy) + t ^ ^ — (15) 

+2{n + A)SV^ = -21) (CyJV^ + CyySVy) + ty^^^ {l6) 



dt Po dx « «/ Pq dx 

,0/0, _ oo/„ ... , „ ... ^ , . c2 dSp 



, , ^^i^'^j; \ ^yx^ '^x , ^yy^ '^y i , ^yx o 

dt Po dx 

+ ik^cl— = -2Vlc^Jv^ + h^^ik^cl— (17) 
dt Po Po , 
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where the quantities with the subscript o are those in the unperturbed stationary disk and is the 
background sound speed. Adiabatic perturbations 

^ = 7^ (18) 
Po Po 

have been assumed, with 7 = 5/3. 

We adopt units of length and time such that i7 = 1 and = Itt. Since we assume a Keplerian 
disk, A = — (3/4)i7. The description of the basic state is completed, apart from the tidal distortion, 
by specifying the sound speed. We have set = 3. Note that as we increase c^kz/^l ~ k^H , inertial 
waves with frequencies ^ i7 become more nearly incompressible. 

The global growth of perturbations depends on the global shape of the tidal distortion. 
Therefore we allow artificial variations in the tidal coefficients (cjj, bi) along the radial (x) direction. 
We consider two different cases: (a) the tidal strength is constant in x and (b) the tidal strength 
has a Gaussian profile, 

- Xo' 



exp 



Aa 



(19) 



with xq = 5 and Axi = 1.5; our computational grids extend from a; = to a; = 10. The strength 
and the oscillation period of the tidal distortion scale by an overall factor that is roughly linear in 
Mg/Mp but has a more complicated dependence on ro/ri,; cf. RG. We fix ro/rb at 0.4 and take 
MjMp in the range 0.05 - 0.2. 

For a given locally homogeneous tide — that is, for (cjj, bi) independent of x — the local growth 
rate is largest for waves whose natural frequencies lie in exact parametric resonance with the tide. 
The meridional wavenumbers (k^jk^) of such modes satisfy 

kl = kl (^-^'-^^ , + ^^-^^^ (20) 

(see RG). For our choices of c^, i7b, and ^l(ro), the corresponding ratio of k^/k^ lies between 0.916 
and 0.955. 

We have used two different radial boundary conditions, refiecting and absorbing. The refiecting 
condition is d^^p = d^^Vy = dx^v^ = Sv^ = 0. The absorbing condition is created by adding a 
term. Since it appears that these waves cannot be absorbed at a single grid point or pair of grid 
points, the absorbing layer is created by adding the term 

, 1 1 a; \ „ . , , 
with - = — 1 - for < a; < Ax2, (21) 

T T To \ AX2J 

to the right hand side of equations (15)-(17). We have used tq = 4 and Ax2 = 2. The round-trip 
attenuation of a wave packet traveling at group velocity Vg from x = Ax2 to a; = and back again 
is 

Aln^t;^-2 / ^^^-5.9, (22) 

J VgT{x) 
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assuming Vg = 0.085 (see below). Note also that kxAx2 ~ 12 for the fastest-growing modes in our 
simulations, so the damping can be estimated with a WKB formula like (22). Thus, the waves of 
interest are almost completely absorbed in the layer < a; < Ax2. 

The linearized perturbation equations (14)-(17) constitute an initial- value problem. We have 
used a localized initial perturbation, 

f Asin [27r(a; - 4.5)1 in 4.5 < a; < 5.5 
[ 0, otherwise, 
Sp = Sp = = Svy = 0, (23) 

with A = 10~^°. We have evolved the equations using the Runge-Kutta-Verner fifth and sixth 
order method with 1000 equally spaced cells along < a; < 10. 



3. Results 

In the calculations with a spatially constant tidal distortion, the perturbations grow exponen- 
tially with localized disturbances slowly spreading out with group velocity (G93) 

1/2 

(24) 

which is ~ 0.085 for our parameters. Fig. 1 shows the temporal evolution of the perturbation 
quantities averaged in quadrature over the grid in a calculation with refiecting boundaries on both 
sides. Here, the strength and the oscillation period of the tidal are appropriate for Mg/Mp = 0.2 
and ro/rb = 0.4, thus i7b/f^ ~ 0.27. From this plot, the measured growth rate is s = 0.139, in good 
agreement with the local growth rate predicted by G93. Fig. 2 shows the velocity eigenfunctions at 
four epochs. Evidently the perturbations spread out not at the sound velocity (c^ = 3) but at the 
much smaller inertial-wave group velocity (24). It is interesting to note that as the perturbations 
grow, the dominant x and z wavenumbers come to have the approximate ratio k^/k^ ~ 0.93, close 
to the one which gives the maximum local growth equation (20). 

When we repeat these calculations with an absorbing layer rather than a refiecting boundary 
on the left side, we see a nearly identical exponential growth rate. In Fig. 3, we have plotted the 
velocity eigenfunctions for a calculation with the same tidal parameters as in Figs. 1&2. The wave 
is clearly suppressed in the absorbing layer, but otherwise it looks very similar to that of Fig. 2. We 
have also done the calculation with absorbing layers at both boundaries (not shown); the effect on 
the growth rate is negligible. Evidently, the instability does not require refiecting radial boundaries. 

Next we take a tidal distortion that varies along the radial direction according to the gaussian 
profile (19). The strength of the tidal terms at the boundaries (a; = and 10) is exp[— (5/1.5)^] = 
1.5 X 10~^ times smaller than the peak strength at a; = 5, and the local growth rate is reduced in 
direct proportion. We take a refiecting boundary on the right and an absorbing layer on the left. 




- 8 - 



From Fig. 4, the global growth rate is s 
than, the rate observed previously. 



0.116, which is comparable to, but significantly smaller 



Fig. 5 shows the velocity perturbations at four different epochs. Since the wave has the same 
shape in the last three epochs, we may consider it a global mode. The eigenfunction is strongly 
peaked at the center of the grid where the local growth rate is largest. However the wings of the 
eigenfunction are clearly growing at the same exponential rate as the central parts, not at their 
own local growth rate. We can fit the eigenfunction to an envelope of the form 



(X exp 



5 



(25) 



with fg ~ 0.8. Notice that this envelope differs in shape from the gaussian profile of the local 
growth rate (19). 

Our interpretation of the results in Figs. 4&5 is that the wave grows locally in the central region 
where the local growth rate is largest; however, because the mode is also a propagating wave, it 
leaks out into the surrounding regions where the local growth rate is negligible. The amplitude at 
time t in the wings Sx from the center is comparable to the central amplitude at the retarded time 
t — \Sx/vg\; this explains the envelope shape (25). 

On this interpretation, the global growth rate should depend on the maximum local growth 
rate, the size of the rapid growth region, and the local group velocity. If disturbances escape from 
the growth region in less than a nominal e-folding time, we expect the instability to be suppressed. 
In order to test and quantify this hypothesis, we have done several calculations with the same 
gaussian profile (19) for the tidal terms, but with different maximum strengths. The latter were 
chosen as appropriate for ro/rb = 0.4 and Mg/Mp = 0.16, 0.13, 0.1, 0.08, and 0.05, as well as 0.2. 
In Fig. 6, we have plotted the global growth rate s, normalized to the maximum local growth rate 
•Smax(?'07 Ms, Mp), agaiust (vg/ Ls^i^^). Here i is a measure of the size of the rapid growth region. 



•Slo 



-dx 



The global growth rate is well fit by 



exp 



1 



Xq 



Axi 



dx ~ 2.66. 



(2Q) 



(27) 



Thus the global growth rate is comparable to the maximum local rate unless Vg ^ Ls^i^^, in which 
case instability is quenched. With appropriate values for L and Smax, this formula may be used 
to estimate the global growth rate in accretion disks. We suspect that in cases where the local 
rate peaks in an extremely narrow annulus, the value of Smax to be used in this formula should be 
obtained by averaging over a radial range Ar ~ l/k^ around the peak. 
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4. Summary and Conclusions 

We have shown that the excitation of perturbations by the elliptical instability in tidally 
distorted disks is essentially local, and is insensitive to boundary conditions as long as the unstable 
region is far enough from the boundaries. If the local growth rate varies radially then regions of 
rapid growth launch waves into the regions of slower growth. The global growth rate is proportional 
to the maximum local growth rate, but is reduced by the ratio of the group velocity to the size of the 
rapid-growth region, i.e. by the rate at which disturbances escape the region. We conclude that the 
tidal excitation of internal waves is insensitive to radial boundary conditions. Since previous work 
has shown that it is only moderately sensitive to vertical structure (Vishniac & Zhang (1995)) it 
follows that this process will inevitably occur in tidally distorted accretion disks. Given the sensitive 
dependence of the growth rate on radius, the torque exerted by the creation of the internal waves 
is likely to play a major role in determining the outer boundary of a disk. 

Our results also imply that it is possible for the m = 2 component of the tidal stresses to excite 
waves with m 1, although not as easily as it will drive waves with m = 1. A wave with some value 
of m, not necessarily equal to 1, will drive a second wave with kr2 = —kri, ^2 = 2(1) — fib) — ^1, 
1712 = 2 — nil, and kz2 = —kzi- The second wave can, in turn, beat against the tidal stresses 
to reinforce the first wave. As long as nii/r and 1112/7 are small compared to k^i and kj-i these 
resonant conditions imply that uji =0)2- Since a) is a different function of r for waves with different 
values of m the resonant conditions can only be satisfied over some annulus of size Ar such that 

ro + Ar 

{kri + kr2)dr < 1. (28) 

If we assume that kj-i + kr2 = in the middle of this annulus then 

dkr » 1 / 3^ (mi — in2)\ , Ar, 

Kx + K2 = -g^Au = ^ (r-ro- — ). (29) 

Equations (28) and (29) imply an upper limit on Ar of 

^r<(4\Y'\ (30, 

where Am = |mi — m2|. Following equation (27) we see that the condition for the existence of an 
elliptical instability involving a wave pair (mi, 2 — mi) is 



SAmfifg 

8r 



This condition can always be satisfied for Vg very small, i.e. for very large wavenumbers; in fact, 
using equation (2) for Smax and (24) for fg, we find when r ^ rb that 

'mA^ f r \ f r \2i/4 



Am < 2kzH — . (32) 
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On the other hand, such waves will have very small saturation amplitudes. If we restrict ourselves 
to large scale modes, i.e. k^H of order unity, then equation (31) is satisfied for Am only in the 
outer parts of thin disks. 

This mechanism drives standing waves, i.e. waves with equal amplitude in outward and inward 
propagating waves. Nevertheless, since the growth rate varies with radius as r"^/^ far from the 
vertical resonance, and much more steeply near it, the local wave field will in general be asymmetric, 
with more waves moving down the local gradient in s than moving up it. This implies that at small 
radii we expect to be dominated by inward propagating internal waves with m = 1, carrying negative 
angular momentum. This will produce a positive angular momentum fiux whose amplitude will 
depend on whether the amplitude of the waves is limited by nonlinear interactions between these 
waves or by interactions with some independently generated local turbulence. Vishniac & Diamond 
1989 assumed the former and estimated an effective (nonlocal) dimensionless viscosity a ~ {H jrY. 
Whether or not this effect is significant depends on the nature of competing processes. 

This work was supported in part by the Non-Directed Research Fund of Korea Research 
Foundation 1993 (DR), and by NASA grants NAGW-2419 and NAG5-2796 (JG) and NAG5-2773 
(ETV). 
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FIGURE CAPTIONS 

Fig. 1 Evolution of spatially averaged r.m.s. velocity perturbations, assuming a spatially constant 
tidal distortion appropriate for Mg/Mp = 0.2, ro/rb = 0.4. Reflecting radial boundary 
conditions. 

Fig. 2 Velocity waveform at four epochs for the same calculation as in Fig. 1. 

Fig. 3 Like Fig. 2, but with a reflecting boundary at the right and and an absorbing layer on the 
left [cf. Eqn. (21)]. Note that amplitude is almost identical with Fig. 2 at corresponding 
times except in the absorbing layer < a; < 2. 

Fig. 4 Evolution of spatially averaged perturbations for tidal distortion with peak strength as in 
Figs. 1-3 but modulated by the gaussian profile (19). Refiecting and absorbing boundaries as 
in Fig. 3. 

Fig. 5 Waveforms at four epochs for the same calculation as in Fig. 4. Note that abscissa is 
logarithmic rather than linear in the final panel; compare equation (25). 

Fig. 6 The global growth rate normalized with the maximum local growth rate (s/smax) as a 
function of the ratio of the group velocity to the size of the rapid growth region times 
the maximum local growth rate (fg/ismax)- Tidal terms modulated by gaussian (19), with 
peak strengths and oscillation periods for Mg/Mp = 0.2, 0.16, 0.13, 0.1, 0.08, and 0.05 and 
ro/rb = 0.4. Boundary conditions as in Figs. 3-5. 
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